********************************************************************************
*Plot all results
********************************************************************************




*plot all survival estimates in a single graph
grc1leg order1_B.gph order2_B.gph order3_B.gph order4_B.gph order5_B.gph, /* 
*/		col(1) subtitle("Original model") saving(col1, replace) /* 
*/		leg(order1_B.gph) ycommon
grc1leg order1.gph order2.gph order3.gph order4.gph order5.gph, col(1) /* 
*/		subtitle("Restricted model") saving(col2, replace) leg(order1.gph) /*
*/		ycommon
grc1leg col1.gph col2.gph, ycommon /* 
*/		xcommon leg(col2.gph) l2("Proportion of cases still at peace") /* 
*/		b2("Days at peace", ring(0) size(small)) ring(5) col(2) 
graph display, ysize(10)
graph export Beardsley_survival_atmean_strata.pdf, as(pdf) replace

